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Abstract 

We search for CP-violating asymmetries (^cp) in the B meson decays to 
K^TT^, K=^7r°, K'^TT^, K'^rj', and ujtt^. Using 9.66 million T(45) decays col- 
lected with the CLEO detector, the statistical precision on ^cp is in the range 
of ±0.12 to ±0.25 depending on decay mode. While CP-violating asymme- 
tries of up to ±0.5 are possible within the Standard Model, the measured 
asymmetries are consistent with zero in all five decay modes studied. 
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CP-violating phenomena arise in the Standard Model because of the single complex parameter 
in the quark mixing matrix [||. Such phenomena are expected to occur widely in B meson decays 
and will be searched for by all current i?-physics initiatives in the world. However, there is currently 
no firm experimental evidence for CP violation outside the neutral kaon system, where direct CP 
violation has been recently observed [|[. 

Direct CP violation, i.e., a difference between the rates for — > / and B ^ f, will occur in 
any decay mode for which there are two or more contributing amplitudes which differ in both weak 
and strong phases. This rate difference gives rise to an asymmetry, Aqp, defined as 

_ B{B ^ /) - B{B ^ /) 

For the simple case of two amplitudes T, P with T <^ P, Acp is given by 

T 

Acp ~ 2| — I sin A(/>wSin A(/)s. (2) 

Here A(ps and A<f)„ refer to the difference in strong and weak phases between T and P. 

The decay B K^tt^, for instance, involves ab u W-emission amplitude (T) with the weak 
phase Arg(V*|_,Vyg) = 7 and a b ^ s penguin amplitude (P) with the weak phase Arg(Vj*'|_,Vjg) = vr 
or Arg(V*bVcJ = @. Theoretical expectations of \T/P\ ~ 1/4 in B ^ K^ir^ g thus allow for 
^CP as large as ±0.5. 

The CP-violating phases may arise from either the Standard Model CKM matrix or from 
new physics while the CP-conserving strong phases may arise from the absorptive part of a 
penguin diagram ^ or from final state interaction effects 0]. Precise predictions for Acp are not 
feasible at present as both the absolute value and the strong interaction phases of the contributing 
amplitudes are not calculable. However, numerical estimates can be made under well-defined model 
assumptions and the dependence on both model parameters and CKM parameters can be probed. 
Recent calculations of CP asymmetries under the assumption of factorization have been published 
by Ali et al. |§ and are listed in Table | for the modes examined in this paper. A notable feature 
of the model used in Ref. Q is that soft final state interactions are neglected, leading to rather 
small CP-invariant phases. However, it has been argued recently that CP-conserving phases due 
to soft rescattering could be large Q, possibly leading to enhanced |^cp| @- 

In this Letter, we present results of searches for CP violation in decays of B mesons to the 
three Ktt modes, K^tt^, K^tt^, Kgir^, the mode K^rj', and the vector-pseudoscalar mode wvr^. 
These decay modes are selected because they have well-measured branching ratios and significant 



signal yields in our data sample [10|. In addition, these decays are self-tagging; the flavor of the 
parent b or b quark is tagged simply by the sign of the high momentum charged hadron. In the 
decay B K^tt^ we assume that the charge of the kaon tags the charge of the b quark. 

The data used in this analysis was collected with two configurations of the CLEO detector at 
the Cornell Electron Storage Ring (CESR). It consists of an integrated luminosity of 9.13 fb^^ 
taken on the T(4S) resonance, corresponding to 9.66 million BB pairs, and 4.35 fb^^ taken below 
BB threshold, used for continuum background studies. CLEO is a general purpose solenoidal 



magnet detector, described in detail elsewhere |11]. Cylindrical drift chambers in a 1.5T solenoidal 



magnetic field measure momenta and specific ionization {dE/dx) of charged tracks. Photons are 
detected using a 7800-crystal CsI(Tl) electromagnetic calorimeter. For the second configuration, 
the innermost tracking chamber was replaced by a 3-layer, double-sided silicon vertex detector, 
and the gas in the main drift chamber was changed from an argon-ethane to a helium-propane 
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mixture. These modifications led to improved dE/dx resolution in the main drift chamber, as well 
as improved momentum resolution. Two thirds of the data used in the present analysis was taken 
with the improved detector configuration. 

Efficient track quality requirements are imposed on charged tracks. Pions and kaons are identi- 
fied by dE/dx. The separation between kaons and pions for typical signal momenta p ~ 2.6 GeV/c 
is 1.7 and 2.0 standard deviations (cr) for the two detector configurations. Candidate are se- 
lected from pairs of tracks forming well-measured displaced vertices with a tt~^tt~ invariant mass 
within 2a of the mass. Pairs of photons with an invariant mass within 2.5a of the nominal 
vr'' mass are kinematically fitted with the mass constrained to the nominal vr*^ mass. For the high 
momentum Kg and vr^ candidates reconstructed with these requirements, the ratio of signal to 
combinatoric background is better than 10. Electrons are rejected based on dE/dx and the ratio 
of the track momentum to the associated shower energy in the Csl calorimeter; muons are rejected 
based on the penetration depth in the instrumented steel flux return. Resonances are reconstructed 
through the decay channels rj' — > r/vr+vr" with rj 77, rj' p7 with p — > 7r"'"7r~, and u — > TT~^ir~Tr^. 

The ^cp analyses presented here are closely related to the corresponding branching fraction 
determinations published elsewhere fin]. We briefly summarize here the main points of the analysis. 



We calculate a beam-constrained B mass M = y E'^ — P^, where is the B candidate mo- 
mentum and -Eb is the beam energy. The resolution in M ranges from 2.5 to 3.0 MeV, where the 
larger resolution corresponds to the B^ K^tt'^ decay. We define AE = Ei + E2 — E\^, where 
El and E2 are the energies of the daughters of the B meson candidate. The resolution on /S.E is 
mode-dependent. For final states without photons, the /S.E resolutions for the two configurations 
of the CLEO detector are 26 and 20 MeV. We accept candidates with M within 5.2 — 5.3 GeV 
and \A.E\ < 200 MeV, and extract yields and asymmetries with an unbinned maximum likelihood 
fit. The fiducial region in M and AE includes the signal region and a substantial sideband for 
background determination. Sideband regions are also included around each of the resonance masses 
{r]' , r], and uj) for use in the likelihood fit. For the ij' pj case, the p mass is not included in the 
fit; we require 0.5 GeV < rriTrTr < 0.9 GeV. 

The main background arises from e~^e~ qq (where q = u,d,s,c). Such events typically 
exhibit a two-jet structure and can produce high momentum back-to-back tracks in the fiducial 
region. To reduce contamination from these events, we calculate the angle Og between the sphericity 
axis [12 1 of the candidate tracks and showers and the sphericity axis of the rest of the event. The 
distribution of cos 6s is strongly peaked at ±1 for qq events and is nearly flat for BB events. For Kir 
modes, we require |cos0s| < 0.8 which eliminates 83% of the background. For r]' and co modes, the 
requirement is |cos0s| < 0.9. detail in Additional discrimination between signal and qq background 
is obtained from event shape information used in a Fisher discriminant {J-") technique as described 
in detail in Ref. [p^] . 

Using a detailed GEANT-based Monte Carlo simulation |l^ we determine overall detection 
efficiencies of 0.48 {K^n^), 0.38 (K^tt^), 0.15 (K°7r=^), 0.13 (K^r]'), and 0.26 (w7r±). These 
efficiencies include secondary branching fractions for Kg vr^vr^ and vr'^ — > 77 as well as 

for the r]' and lo decay modes where applicable. 

To extract signal and background yields we perform unbinned maximum-likelihood fits using 
AE, M, [cos (if not used in T\ dE/dx, daughter resonance mass, and helicity angle in 
the daughter decay. The free parameters to be fitted are the asymmetry (/ — /)/(/ + /) and the 
sum (/ -|- /) in both signal and background. In most cases there is more than one possible signal 
hypothesis and its corresponding background hypothesis, e.g., we fit simultaneously for K^t^^ and 
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IT TT^ to ensure proper handling of the K/tt identification information. The probabihty density 
functions (PDFs) describing the distribution of events in each variable are parametrized by simple 
forms (Gaussian, polynomial, etc.) whose parameters are determined in separate studies. For 
signal PDF shapes parameters are determined by fitting simulated signal events. Backgrounds 
in these analyses are dominated by continuum e~^e~ — > qq events, and we determine parameters 
of the background PDFs by fitting data collected below the T{4S) resonance. The uncertainties 
associated with such fits are charge symmetric in all PDFs except the dE/dx parametrization. The 
dE/dx information was calibrated such that any residual charge asymmetry is negligible compared 
to the statistical errors for Ac^. 

The experimental determination of charge asymmetries in this analysis depends entirely on the 
properties of high momentum tracks. The charged meson that tags the parent h/h flavor has in each 
case a momentum between 2.3 and 2.8 GeV/c. In independent studies, using very large samples of 
high momentum tracks, we searched for and set stringent limits on the extent of possible charge- 
correlated bias in the CLEO detector and analysis chain for tracks in the 2 — 3 GeV/c momentum 
range. Based on a sample of 8 million tracks, we find an ^cp bias of less than ±0.002 introduced 
by differences in reconstruction efficiencies for positive and negative high momentum tracks. 

For K^TT^ combinations, where differential charge-correlated efficiencies must also be consid- 
ered in correlation with K/-K flavor, we use 37,000 — > i^7r(7r'') decays and set a limit on the ^cp 
bias of ±0.005. These D'^ meson decays, together with an additional 24,000 Z?^^ meson decays, 
are also used to set an upper limit of 0.4 MeV/c on any charge-correlated or charge-strangeness- 
correlated bias in the momentum measurement. The resulting limit on Acp bias from this source 
is ±0.002. We conclude that there is no significant ^cp bias introduced by track reconstruction 
or selection. 

Our ability to distinguish the final states K^tt^ and K^tt^ depends entirely on particle iden- 
tification using dE/dx. In addition, all other decay modes depend to varying degrees on dE/dx to 
distinguish between B Xtt^ and XK^, X being a Kg,7r^,ri' or an lo. The dE/dx was carefully 
calibrated in order to remove any possible charge dependencies. 

We calibrate the dE/dx response using radiative // pair events assuming that for a given 
velocity dE/dx is the same for /i^, vr^, and K^. We then compare the dE/dx response for 
positive and negative tracks in the momentum range 2 — 3 GeV/c from all hadronic events in the 
CLEO data sample. The large available statistics allows us to split the data into subsets and to 
verify the stability of the calibration over time. The fully calibrated dE/dx is then verified using 
kinematically identified kaons and pions of 2 - 3 GeV/c from K'k{'k^) decays. The dE/dx 

distributions for and vr^ from this sample are shown in Fig. 0. No significant differences are 
seen between different charge species. The statistical uncertainty in this comparison translates into 
a possible ^cp bias of ±0.01 for K^tt"^ , and less for all other final states. We conservatively assign 
a total systematic error of ±0.02 in all five decay modes. 

As additional check we measure the asymmetry of the background events in each decay mode, 
and find that all are consistent with the expected null result for continuum background. The 
results for the asymmetry in continuum background are -0.024 ±0.038 (K=^7r^), -0.003 ± 0.032 
(i^±7r0), -0.017 ± 0.037 (i^gvr±), -0.006 ± 0.070 (V(??vr7r)K±), -0.009 ± 0.015 {rj'{pj)K^), and 
-0.001 ± 0.010 {ujTT^). We further confirm that our analysis method does not introduce a bias in 
the measured ^cp in the analysis of simulated events with known asymmetries. 

We conclude that any possible systematic bias on ^cp is negligible compared to the statistical 
errors of our measurements. Our 90% confidence level (CL) ranges are calculated adding statistical 
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and systematic errors in quadrature. 



TABLE I. Summary of results. Signal yields are taken from Ref. 10. Theory predictions are 
from Ref. 8, and include only Standard Model perturbative calculations. The 90% CL interval 
includes statistical and systematic errors (±0.02) added in quadrature. 



Mode 


Signal 


^CP 


^CP 


^CP 




Yield 




90% CL 


Theory 




80+J^ 


-0.04 ±0.16 


[-0.30,0.22] 


(±0.037, ±0.106) 




42.lt^V 


-0.29 ±0.23 


[-0.67,0.09] 


(±0.026, ±0.092) 




ot; 9+6.4 


+0.18 ±0.24 


[-0.22,0.56] 


±0.015 




lootii 


±0.03 ±0.12 


[-0.17,0.23] 


(±0.020, ±0.061) 




OQ c;+8.2 


-0.34 ±0.25 


[-0.75,0.07] 


(-0.120, ±0.024) 



We summarize the results in Table | and Fig. The dependence of the likelihood function 
on ^cp for each of the five decay modes is depicted in Fig. ^ This figure was obtained by re- 
optimizing the likelihood function at each fixed value of ^cp to account for correlations between 
the free parameters in the fit. 

We see no evidence for CP violation in the five modes analyzed here and set 90% CL intervals 
that reduce the possible range of ^cp by as much as a factor of four. It has been suggested |15] 
that ^cp in K^ir^ and K^ir^ are expected to be almost identical within the Standard Model. 
Based on the average ^cp in these two decay modes we calculate a 90% CL range of —0.28 < ^cp 
< ±0.05. 

While the sensitivity is not yet sufficient to probe the rather small ^cp values predicted by 
factorization models, extremely large Acf values that might arise if large strong phase differences 
were available from final state interactions are firmly ruled out. For the cases of Ktt and r]'K, we 
can exclude |^cp| greater than 0.30 and 0.23 at 90% CL respectively. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent lumi- 
nosity and running conditions. This work was supported by the National Science Foundation, 
the U.S. Department of Energy, the Research Corporation, the Natural Sciences and Engineering 
Research Council of Canada, the A.P.Sloan Foundation, the Swiss National Science Foundation, 
and Alexander von Humboldt Stiftung. 
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FIG. 1. Comparison of normalized dE/dx distributions for and (above); and vr^ and 
7r~ (below). The inset shows the asymmetry {K~ — K^)/{K~ + K^) and (vr" — 7r"'")/(7r~ + vr"'") 
respectively. Kinematically identified kaons and pions of momenta 2 — 3 GeV/c from D° K'k{'k^) 
decays in data are used for this comparison. 
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FIG. 2. measurements. Error bars include systematics; hatched regions delimit the yu7o 

CL intervals. 
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